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Transmission Capacity of Wireless Ad Hoc 
Networks with Energy Harvesting Nodes 



Rahul Vaze 



Abstract 

Transmission capacity of an ad hoc wireless network is analyzed when each node of the network 
^ . harvests energy from nature, e.g. solar, wind, vibration etc. Transmission capacity is the maximum 

allowable density of nodes, satisfying a per transmitter-receiver rate, and an outage probability constraint. 
Energy arrivals at each node are assumed to follow a Bernoulli distribution, and each node stores energy 
using an energy buffer/battery. For ALOHA medium access protocol (MAP), optimal transmission 

in 

probability that maximizes the transmission capacity is derived as a function of the energy arrival 

o 

. distribution. Game theoretic analysis is also presented for ALOHA MAP, where each transmitter tries 

to maximize its own throughput, and symmetric Nash equilibrium is derived. For CSMA MAP, back-off 
probability and outage probability are derived in terms of input energy distribution, thereby characterizing 

&• 

the transmission capacity. 



I. Introduction 

Consider an ad hoc wireless network, where multiple source- destination pairs try to commu- 
nicate without the help of a centralized controller. There are many examples of real life ad hoc 
networks, such as military networks, vehicular networks, sensor networks, etc. Typically, each 
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node of an ad hoc network is powered by a conventional battery, that has limited lifetime and 
needs to be replenished periodically. Recently, to improve the lifetime of nodes, and to provide a 
means of green communication, the concept of equipping nodes with energy harvesting devices 
that can extract or tap energy from renewable energy sources, such as solar, wind, vibration, 
etc., has been introduced. Harvesting energy from nature, however, makes the future energy 
availability random, and the transmission strategies have to adapt dynamically to the energy 
arrivals. 

In this paper, we are interested in characterizing the transmission capacity of an ad hoc network 
when each transmitter has energy harvesting capability. The transmission capacity of an ad hoc 
network is the maximum allowable density of nodes, satisfying a per transmitter-receiver rate, 
and outage probability constraint jT), ED. Similar to (Tj, Q, we assume that the locations of 
transmitters are distributed as a homogenous Poisson point process (PPP) with a fixed density. We 
assume that energy arriving to each node through natural sources follows a Bernoulli distribution, 
where in each discrete time, a unit amount of energy arrives with probability p or no energy 
arrives with probability 1 — p. Energy arrivals are independent and identically distributed across 
all transmitters in the network. We remark in the sequel that even though we consider Bernoulli 
energy arrival model, our results apply to general energy arrival distributions with unbounded 
battery capacity. 

In this paper we consider two medium access protocols (MAPs): ALOHA, and CSMA, to be 
used by each transmitter. For ALOHA MAP, we derive the optimal transmission probability as 
a function of energy harvesting rate p that maximizes the transmission capacity of the ad hoc 
network. We also consider a selfish setting, where each transmitter is interested in maximizing 
its own throughput, and derive symmetric Nash equilibrium strategies. For the CSMA MAP, 
we derive the back-off probability and the outage probability that together characterize the 
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transmission capacity. Before describing our contributions, we first review some recent work on 
the design of wireless networks with energy harvesting. 

A. Prior Work 

Assuming non-causal knowledge of future energy arrivals, optimal offline power allocation 
strategies to maximize throughput have been derived for single source-destination pair [|3), 
interference channel [4], and broadcast channel [5). For a single source-destination pair, structure 
of a causal throughput optimal strategy with energy arrival distribution knowledge has been 
derived in [|6l, while a distribution free online algorithm has been derived in ||7). For a wireless 
network with energy harvesting nodes, the probability that a node successfully transmits a data 
packet to a fusion center is analyzed for time division multiple access and ALOHA MAP in 
lOO. From a queuing theoretic point of view, queue stabilizing policies have been derived for a 
single source-destination pair |]9), and a two-user interference channel IfTOl . 

The most relevant work to this paper is [fTD . where each node of the ad hoc network harvests 
energy with arbitrary energy arrival distribution with fixed rate of arrival p. Each transmitter is 
scheduled to transmit with power P if it has more than P amount of energy irrespective of all 
other nodes. Optimal value of P is derived in [fTTI that maximizes the transmission capacity as 
a function of network parameters. Compared to [11], in this paper we take a different viewpoint 
and couple the energy queue dynamics with the ALOHA/CSMA transmission probability, and 
try to find the best ALOHA transmission probability that maximizes the transmission capacity. 
In our model, a transmitter sends a packet with probability q with ALOHA (if channel is sensed 
idle with CSMA) if there is energy available at the transmitter. 

For the selfish setting, where each transmitter of an ad hoc network powered with a con- 
ventional power source tries to unilaterally maximize its own throughput, optimal transmission 
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probability for ALOHA MAP has been derived in lfT2l . Not surprisingly an always transmit 
strategy is selfishly optimal, and to improve the selfish behavior towards globally optimal 
solution, a penalty function is introduced in the objective function that linearly increases the 
penalty with the transmission probability of each node. 

B. Contributions 

• For ALOHA MAP, we find the optimal transmission probability that each transmitter should 
use to maximize the transmission capacity of the ad hoc network for both finite and infinite 
battery capacity. For the infinite battery capacity case, the optimal transmission probability 
takes two values depending on a threshold that is a function of the network parameters. 
For finite battery capacity case, the optimal transmission probability is shown to satisfy the 
probability of having non-zero energy in the energy queue to be equal to a constant, with 
explicit solution found for the unit battery capacity case. 

• For ALOHA MAP, we find the optimal transmission probability that each transmitter should 
use that maximizes its own throughput (selfish setting), and characterize symmetric Nash 
equilibrium (SNE). At SNE, each transmitter uses transmission probability q such that 
V < Q < 1 w i m infinite battery capacity, and q = 1 with finite battery capacity. It is 
shown that the price of anarchy, that is ratio of globally optimal transmission capacity to 
transmission capacity obtained at the worst SNE, is quite low and is actually equal to one for 
some cases. Thus, in contrast to |[T2l . the selfishly optimal strategy in the energy harvesting 
setting is not much different from the globally optimal strategy, since with energy harvesting 
nodes, individual objective functions are inherently energy aware. 

• For CSMA MAP, we derive the back-off probability and the outage probability for each 
transmitter when each transmitter harvests energy from renewable sources, consequently the 
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transmission capacity. 

We remark that using a simple Bernoulli energy arrival model and energy queue based 
transmission probability framework, this work takes few initial steps towards understanding the 
fundamental performance of energy harvesting ad hoc networks, and a lot more work is required 
for complete characterization. 

II. System Model for ALOHA MAP 

Consider a wireless ad hoc network where multiple source destinations pairs want to commu- 
nicate with each other without any centralized control. Following HI, the location of transmitter 
nodes T m , m E N is assumed to be distributed as a homogenous Poisson point process (PPP) 
$ = {T m } on a two-dimensional plane with density A [fT3l . The receiver R m associated with 
transmitter T m is assumed to be at a fixed distance of d from T m with an arbitrary orientation. 

We assume that each transmitter harvests energy from nature, e.g. through solar, wind, peizo- 
eletric sources etc. Each transmitter is assumed to have a battery of capacity B using which it 
can store the harvested energy. The energy arrival process is assumed to be i.i.d. Bernoulli with 
rate p across different transmitters, i.e. at each time t, either a unit amount of energy arrives at 
any transmitter with probability p, or no energy arrives with probability 1 — p. The results of 
this paper are applicable to general energy arrivals as pointed out in Remark Q] The Bernoulli 
energy arrival assumption is made for simplicity of exposition and for deriving critical insights 
into the problem. 

If x m is the transmitted signal from T m at time t, then the received signal y m at R m is given 

by 

y t m = VPd- a / 2 h mm x m + v / Pl Ta (t)rf~: /2 ^ s + ^ m , (1) 

r s e*\{T } 



where P is the power transmitted by each transmitter, a > 2 is the path-loss exponent, l Ts (t) 
is the indicator function that is 1 if T s transmits at time t and is zero otherwise, hu is the 
fading channel coefficient between transmitter Ti and receiver Rk that is assumed to be Rayleigh 
distributed, and z m is the zero mean unit variance additive white Gaussian noise. 

We consider the interference limited regime, i.e. noise power is negligible compared to the 
interference power, and henceforth drop the noise contribution [Q]|tj With the interference limited 
regime, we assume unit power transmission, P — 1, since the signal to interference ratio (SIR) 
does not depend on P. Let SIR m (t) denote the SIR between T m and R m at time t, then using 
© 



SIRm(t) : = 



d~ a \h I 2 



12 



lT s (t)d- a s \h r 

We consider a slotted ALOHA like random access MAR where each transmitter attempts 
to transmit its packet with an access probability q, independently of all other transmitters if it 
has energy to transmit. Thus, with transmit power P = 1, any transmitter attempts to transmit 
if it has non-zero amount of energy. Let E m (t) be the amount of energy available with T m 
at time t. Then E m (t) is an i.i.d. birth-death Markov process with the following transition 
probabilities, P(E m (t + 1) = E m (t) + 1) = p{l - q), P(E m (t + 1) = E m (t) - 1) = q(l - p), 
P{E m {t + 1) = l\E m (t) = 0) = p and P{E m {t + 1) = 0\E m (t) = 0) = 1 -p as shown in Fig. 
ED Let r := P(E m (t) > 1) be the probability of E m (t) > 1. Let Xm(t) = 1 if T m transmits at 
time t given that E m (t) > 1, and Xm(t) = otherwise. Thus, P(xm(t) = l\E m (t) > 1) = q. 
By definition, any transmitter transmits (1t s (*) = 1) with probability P(E m (t) > l,Xm — 1) — 
P(E m (t) > l)P(xm(t) = l\E m (t) > 1) = P(E m (t) > l)q = rq, independently of all other 
nodes. Consequently, the active transmitter process is a homogenous PPP on a two-dimensional 

'This assumption is made for simplicity of exposition, and all results of this paper can be easily extended to the case of 
additive noise as well. 
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plane with density A a := qrX. 

We define that the transmission from T m to R m at time t is successful if the SIR between T m 
and R m is greater than a threshold 8, i.e. SIR m (£) > 8. Thus, the probability of success at time 
t is defined to be 

P suc (t) = P(S\R m (t) > 8), (2) 

We assume that the rate of transmission corresponding to threshold 8 is R = log(l + 8) 
bits/sec/Hz. Then the transmission capacity \X\ is defined to be C = \ a P suc (t)R bits/sec/Hz/m 2 . 

Our goal is find the optimal q that maximizes the transmission capacity, q* = argmax 9 C. For 
the purpose of analyzing the success probability P suc and C, we consider a typical transmitter 
receiver pair T m , R m . It has been shown in [Q]| that for the PPP distributed transmitter locations, 
the performance of the typical source destination pair is identical to the network wide perfor- 
mance using Slivnyak's Theorem |fl3l . Next, we first analyze the case when each transmitter has 
an unbounded battery capacity B = oo, and later extend it to the finite battery capacity case. 

A. B = oo (Infinite Battery Capacity) 

From Fig. [TJ using well-known analysis of the infinite state birth-death Markov process we 
have that r = P(E m (t) > 1) = min j |, 1 j for p, q > 0, V m. Throughout this paper, the only 
quantity we will be interested in from the energy queue point of view is the probability that the 
energy queue is not in state 0, P(E m (t) > 1). Next, we remark that because of this restricted 
dependence, the results of this paper generalize to any energy arrival distribution, and are not 
restricted to just Bernoulli distribution. 

Remark 1: Let the energy arrival process be A(t) with an arbitrary distribution of rate p that 
is independent and identically distributed for all transmitters. Let each transmitter transmit with 
probability q, if there is more than P amount of energy in the battery, and uses power P for its 



transmission (similar to the ALOHA setup described before). The energy queue dynamics at any 
transmitter is given by E(t + 1) = E{t) + A(t) - Pl E{t) > P x(t), where q = E{x(t)\E(t) > P}. 
For this general model with B = oo, it has been shown in |fTTfl that P(E{t) > P) = 1 if p > Pq, 
and P(E(t) >P) — -p-ifp< Pq. In ifTTI . the result is obtained for q — 1, but it can be readily 
generalized for any q e [0, 1]. Specializing, this result for our model with unit power P = 1, 
we get that P(E(t) > 1) = min |~, 1 j for any energy arrival distribution. Thus, we do not lose 
any generality by restricting ourselves to the Bernoulli energy arrival distribution. 
Consider the success probability ©, 

P suc = P(S\R m {t)>6), 

d~ a \h I 2 

p I " V L mm\ q 



g$\ { To} 1 T.(*)dmSl>l 
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exp (— \rqd 2 9 ^k{q)\ , lfl4ll where n(a) 



2tt 2 



asin{2ir/a) ' 



Let Xmax '■= d2s l - - • Hence the transmission capacity is C = Xrqexp ^— -f^-j R- We derive 
the optimal transmission probability q* in Theorem CD We need the following definition for the 
proof of Theorem [T] 

Definition 1: A function / : R — > R is called unimodal if for some value m, it is monotonically 
increasing for x < m and monotonically decreasing for x > m. 

Theorem 1: The optimal ALOHA transmission probability with B = oo that maximizes the 
transmission capacity is q* = if p > ^aa, and any q E [p, 1] is optimal otherwise. 
Proof: With r = min ||, 1 j, for | < 1, the transmission capacity is C = Apexp ^— \ R, 
while with 2 > 1, where r = 1, C = Agexp (-j^J 

Case 1: < 1 and p > ^p.} Note that ^f^ = arg max, \q exp ("a 2 ^ - ) ^ More- 

over, since < 1 andp > ^p, considering < q < p, where r = 1, = argmax < g < p C - 
arg max <g< p g exp ^— y^-j R, since g = lies in the feasible set < q < p. With 



f = , the optimal transmission capacity is C = 



p+q-pq ' 

B 



while for B > 1, r B = v p(p[1 - q) ,B J for p, q > 0,p ^ q and r B = 1TT — for p = q. Similar 



Case 2: Let p < ^saa. Since C = Agexp ^— i? is a unimodal function of q, and 
achieves its maxima at it implies that Agexp ( — l^ - ) ^ ls an increasing function of q for 
< ? < p for p < ^^l. Hence max < 9 < p C = Apexp ^— -^—j R- Moreover, for any q > p, 
C = \pexp ^— R- Hence, if p < ^f^-, then any q £ [p, 1] is optimal. ■ 

Next, we consider the more realistic case of finite battery capacity B at each node. 

B. Finite B (Finite Battery Capacity) 

With finite battery capacity B, the energy queue transition probabilities are illustrated in Fig. 
|2l and let r B := P(E m (t) > 1). With finite battery capacity B, the energy queue is a finite state 
birth-death Markov process for which we have that for B — 1, T\ — P(E m (t) > 1) 

p( i ( p(l-q) \ B \ 

q y L \g(l-p)) J 

to B = oo case, the transmission capacity expression for finite B is C = \r B qexp ^ x , . 
and we want to maximize C with respect to q. We first prove this intermediate result that is 
important for subsequent analysis. 

Lemma 1: For finite B, function f B (q) '■= qr B is an increasing function of q for q £ [0, 1]. 
Proof: See Appendix [A] ■ 

Theorem 2: The optimal ALOHA transmission probability with finite battery capacity B is 
q* = min {q, 1}, where q is the solution to the equation f B (x) = ^aasi. 

Proof: From Lemma [TJ we know that f B (q) is an increasing function of q for q £ [0, 1]. 
Hence for fixed A and X max , C = A f B (q) exp ( — y^— f B (q) ) is a unimodal function of q. 
Thus, either the maxima of C lies in [0, 1], or C is an increasing function for [0, 1]. To find 
the maxima of C, we equate the first derivative of C, C = \f' B (q)exp ^— x^ - /b(<?)) — 
A M<?)A^(<?)exp (-^^(g)) to zero, which yields = 1 and f B (q) = Let q 
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solve }b{q) — ^f^. If q > 1, then q* = 1, otherwise q* = q. ■ 

Corollary 1: With unit battery capacity 5 = 1, q* — min j Ap+ ^ Ama ^_ p p 1 j. 
Proof: From Lemma [j] for S = 1, /{(g) > for g G [0,1]. Hence using Theorem [2l we know 
that q* satisfies A { l(q * } = 1, where A(g) = —22—. Consequently q* = p a ^* — B 

Discussion: In this section, we derived the optimal transmission probability for ALOHA MAP 
when each transmitter is harvesting energy with a rate p Bernoulli distribution. We also remarked 
that our results apply to any general energy arrival distribution. For the infinite battery capacity 
case, we showed that the optimal transmission probability takes two values depending on a 
threshold that is function of the network parameters. There is a natural connection between the 
optimal transmission probability with nodes that are powered with energy harvesting sources 
and conventional power sources [0. If rate p is more than the optimal transmission probability 
with conventional power sources, then we know that the probability that there is more than unit 
energy in the queue is 1. Hence there is always energy to transmit, and the optimal transmission 
probability with energy harvesting case is equal to the conventional power sources. If rate p is 
less than the optimal transmission probability with conventional power sources, then we show 
that transmission capacity is an increasing function of transmission probability from to p 
and then becomes a constant for any transmission probability greater or equal to p. Thus, any 
transmission probability greater or equal to p is shown to be optimal. 

We also derived results for the realistic case of finite B. We showed that the transmission 
capacity expression is unimodal, and its maxima as a function of transmission probability can 
be found by solving the probability of non-zero energy in queue to be equal to a constant. For 
the special case of B = 1, we derived explicit solution for the optimal transmission probability. 

In this section, we considered maximizing system wide transmission capacity that captures the 
sum throughput of the network. In the next section, we consider the case when each transmitter 
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selfishly tries to maximize its own throughput and derive symmetric Nash equilibrium strategies. 



In this section, we consider the case when each transmitter wants to maximize his own 
throughput, and derive the selfishly optimal transmission probability for each transmitter us- 
ing ALOHA MAP. The energy harvesting and transmission protocol at each transmission is 
assumed identical to the previous section except for the transmission probability q, which is 
now transmitter dependent, and we denote the transmission probability of transmitter T n by q n . 
Therefore the objective function (utility) that each transmitter maximizes is it own throughput 
TH n := r n q n P(S\R n (t) > 6)R, where r n = P(E n (t) > 1). 

In this selfish setting, we will consider symmetric Nash equilibrium (SNE), at which all nodes 
use the same transmission probability q*, since no two transmitters are distinguishable from each 
other. For more details see [fT5l . In this setting, let the n th transmitter use q n = q (r n = r), while 
the m th transmitter uses q m = q,m ^ n (r m = f). Then we denote the throughput of the 
n th transmitter as TH n (g, q) := rqP(S\R n (t) > 9)R. Then q* is a SNE if for each transmitter 
TH(g*,g*) = max 9g [ 0i i] TH(g, q*). The transmission capacity of the ad hoc network is defined 
as the sum of throughput of all nodes C = XJH(q*,q*)R. 

A. Infinite Battery Capacity B = oo 

Theorem 3: In the infinite battery capacity case (B = oo), any q* such that p < q* < 1 is 
a SNE, where p is the rate of energy arrivals. Moreover, at any SNE, each transmitter gets the 
same throughput TH = pexp ( — ), and the transmission capacity is C = Xpexp ( —^ L - ). 



III. Selfishly Optimal Transmission Strategy 
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Proof: Consider 

TH n {q,q) = rqP(S\R n (t)> 



d~ a \h I 2 

rqP | ^ > I ' where 1t *W = 1 W "P- ^' 



1 ?r/s 



rgexp 



A 



max 



We note that the throughput of the n th transmitter is monotone increasing in rq. Hence we need 
to find the optimal q that maximizes rq. From the previous section, we know that for the infinite 
battery capacity case {B = oo), r = min ||, 1 1. Thus, if q < p, then r = 1 and rq = q, while 
if q > p, r — - and rq = p. Thus, any q such that p < q < 1, maximizes rq, and provides the 
selfishly optimal throughput for the n th transmitter. Therefore, each q* such that p < q* < 1 is 
a SNE, and the throughput of each transmitter for p < q* < 1 is TH = nexp ( — — ) • ■ 
Since any q such that p < q < 1 is a SNE, it might appear that some transmitters can use q — 1, 
and create more interference for other transmitters. Each transmitter, however, gets to transmit 
only with probability rq — p for p < q < 1, since if q is large, transmitter uses up more energy 
and there is higher chance of energy queue to be in state 0. Next, we compute the price of 
anarchy that compares the performance loss incurred due to selfishness by each transmitter. 

Definition 2: The price of anarchy (PoA) of a game is the ratio of the utility at the globally 
optimal solution to the utility at the worst equilibrium. 

Lemma 2: The PoA of the throughput game is 1 if p < and Am /"i pA n , otherwise. 
Proof: If p < ^n*aa, the globally optimal transmission probability q* is that such that p < q* < 1 
(Theorem [B, which is identical to the selfishly optimal transmission probability q* achieving 
the SNE (Theorem |3]). Hence PoA is 1 when p < ^f^- For V > ^f^' tne globally optimal 
transmission probability is q* = ^aa, whereas the selfishly optimal policy is p < q* < 1. The 
required expression is obtained by taking the ratio of globally optimal transmission capacity 
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expression and transmission capacity expression at any SNE. ■ 
Remark 2: Recently, selfishly optimal transmission probability for ALOHA MAP has been 
derived in [|T2l for throughput maximization (TH(g, q) defined above) with PPP distributed 
transmitters, where each transmitter is powered by a conventional energy source. As expected, 
without any energy constraints, q — 1 (always transmit strategy) is selfishly optimal, however, it 
provides a very poor PoA performance. To improve the selfish behavior towards globally optimal 
solution, a modified objective function is considered TH(g, q) — pq in lfT2l that linearly penalizes 
the increase in transmission probability q, where p is a constant. It has been shown that by 
carefully choosing the scaling parameter p, the PoA can be significantly improved. 

In comparison to [12], for the setting in this paper, where each transmitter harvests energy from 
nature and stores it in a battery, the objective function TH(g, q) is already energy aware, and no 
extra energy dependent factors need to be introduced for obtaining good PoA performance. More 
importantly, we notice that PoA obtained with the improved energy penalty strategy [fl2l is worse 
than the PoA obtained with our model of energy queue dependent transmission probability. Thus, 
even for conventioanal energy powered sources, using a virtual energy queue based transmission 
probability can improve the PoA performance. 

B. Finite B 

Theorem 4: With finite battery capacity B, q* = 1 is a SNE. At any SNE, each transmitter gets 
the same throughput TH = pexp ( — ) and transmission capacity is C = Xpexp ( —t^— J . 
Proof: Similar to the B = oo case, we have that TH n (g, q) = rqexp (^—rqexp (x^ - ))' an( * 
we need to find the optimal q that maximizes rq, i.e. find q that maximizes rq = fsiq)- The 
function fsiq) is an increasing function of q (Lemma [B, and achieves its optimal value equal 
to p at q = 1. Thus q* = 1 is a SNE with throughput of each transmitter TH = pexp ( — y 2 ^- 
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Next, we compute the PoA for the finite battery capacity case. 

Lemma 3: The PoA of the throughput game with finite battery capacity is 1 if y* > 1, where 

. \q* cxp ( — , g A ) . 

y* is the solution to fsiy) = x^, anc * r ™*\ where /s(g*) = otherwise. 

Proof: If satisfies = ^f 21 and y > 1, then the globally optimal transmission probability 

q* = 1 (Theorem |2l) is equal to the selfishly optimal transmission probability q* = 1 (Theorem H]). 
Otherwise, the globally optimal transmission capacity is \q* exp ^— jf^-J , where fsio*) = ^f- 
(Theorem |2]), while the transmission capacity at SNE is Apexp ( — -r 2 ^- ) (Theorem H). ■ 

Discussion: In this section, we considered the game theoretic setting for ALOHA MAP 
where each transmitter unilaterally tries to maximize its own throughput. With each transmitter 
harvesting energy from nature at a finite rate, we showed that the selfishly optimal and globally 
optimal strategies are not very different, and the PoA is quite small. With conventional energy 
powered transmitters, the selfish strategy is to always transmit, however, in the energy harvesting 
setting there is no incentive for any transmitter to transmit aggressively, since more transmission 
attempts deplete the energy available for future transmission. 

In the previous two sections we analyzed the transmission capacity of an ad hoc network with 
ALOHA MAP. Another widely used MAP is CSMA, and in the next section we analyze the 
performance of CSMA in an ad hoc network with energy harvesting nodes. 

IV. CSMA 

In this section, we consider the CSMA MAP and consider a slightly different network model 
compared to Section HH that has been introduced in lfT6l . [fTTl . In Section HH we assumed that 
transmitter locations are distributed as a 2-D PPP, and each transmitter always had a packet to 
transmit and uses ALOHA MAP for packet transmissions. With CSMA, analyzing such model is 
rather challenging. In this section, we consider an area A, and model the packet arrival process 
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as a one-dimensional PPP with arrival rate (A/L)X, where L is the fixed packet duration. Each 
packet after arrival is assigned to a transmitter location that is uniformly distributed in area 
A, and the receiver corresponding to a particular transmitter is located at a fixed distance d 
away with a random orientation, as shown in Fig. |3] For A — > oo, with ALOHA MAP, this 
process corresponds to a 2-D PPP of transmitter locations with density A (Section HO), where 
each transmitter has packet arrival rate of j. 

The SIR between transmitter T n and its receiver R n at time t is SIR n (t) = = — \ hn "}. a — ^, 

Z J r s e*\{T„} 1t s \nna\ 

similar to Section HH where $ is the set of all transmitters, and 1t s (£) = 1, if the transmitter T s 
is not in back-off and has energy to transmit, and otherwise. With CSMA MAP, transmitter T n 
sends its packet at time t if the channel is sensed idle at time t, which in our case corresponds 
to SI R ri (t) > 6, with unit power if available energy E n {t) > 1. Otherwise, the transmitter backs 
off and makes a retransmission attempt after a random amount of time. If T n transmits the 
packet, the packet transmission can still fail if SIR„ falls below 9 for the duration of packet 
transmission L. Thus, the outage probability P out = P b + (1 — Pfc)Pf a ii|no backoff, where P b is the 
back off probability, and Pfaiibo backoff is me probability that the transmission fails. Hence, the 
transmission capacity with CSMA MAP is defined as C = A(l — P ou t)R bits/sec/Hz/m 2 . 

Similar to Section HH we assume that the energy arrival process is i.i.d. Bernoulli with rate 
p across different transmitters. In this section, we only consider the B = oo case. Analysis for 
finite B follows similarly. The transition probability diagram for energy queue with CSMA is 
identical to Fig. [TJ with q replaced by 1 — P&, and r = P(E n (t) > 1) = min |x^f£' 

Remark 3: CSMA MAP introduces correlation among different transmitter's back-off events, 
and hence the number of simultaneously active transmitters on the 2-D plane no longer follows a 
PPP. Nevertheless, for analytical tractability, as an approximation we assume that the transmitter 
back-off events are independent, and simultaneously active transmitter locations are still PPP 
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I, 



distributed. The simulation results show that this assumption is reasonable [fl6l , IfTTI . 

In the next Theorem we derive the back-off probability for any transmitter with the CSMA MAP. 

Theorem 5: The backoff probability P b = 1 - exp (-3^) if ~ Xm ^ lnp > p, otherwise P? 
satisfies P b = 1 — exp ( — x ^~ Pb > ) which can be solved using Lambert's function Wo(-)- 
Proof: Transmitter T n goes into backoff at time if SI R n at time is less than 9. Note that 
the set of transmitters active at time are those that started transmitting between — L to since 
the packet length is L. The transmitters that become active at any time t between time — L and 
is a PPP with density ^(1 — P b )r. Assuming independent back off events across different 
transmitters, the active set of transmitters at any time slot between — L and are independent, 
and since the union of independent PPPs is also a PPP with sum of the densities, transmitters 
that are active at time is a PPP with density Y^,=-l tO- ~ ^ r = ^0- ~ -^>) r - F° r l^ge A this 
translates to having PPP distributed active transmitter locations on the 2-D plane with density 
A(l - P b )r. Thus, P b = P(SIR n (0) < 9) = 1 - exp (-^^ L ) 0. 

Next, we proceed using contradiction. Let P b > 1 — p. Then r = 1, and hence P b = 1 — 
exp ( — ^y— — ) > 1 — p, which results in 1 — P b > ~ Xmax lnp . However, p, A, and A max are fixed 
parameters and if they satisfy the relation ~ Xma ^ lnp > p t it implies that P b < 1 — p. Thus, we 
get a contradiction, since we started with P b > 1 — p. Hence if ~ Xma ^ lnp > p, P b < 1 — p, and 
correspondingly r = and P b = 1 — exp ^— J. The other case is obvious. ■ 



Next, we derive an explicit expression for packet failure probability with the CSMA MAP. 

, l + d a 9x 



|no backoff — 1 : : jTTn • F° r L — 1, 



Theorem 6: P fail | no backoff = 1 

PfaUIno backoff = 1 - (1 - P 6 ) exp (2\9 2 / a d 2 (l - P b )W^CSC (f )). 

Proof: Note that Pf a ii|no backoff is the probability that at any time t, SIR n (£) < 9 for < t < L 
given that SIR n (0) > 9. The transmitters that become active at any time t between time and 
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L is a PPP with density £(1 - P b )r. Then, 




1 



P(SIR n (l) > 9, . . . , S\R n (L) > 9\S\R > 6) 



1 



P(SIR > 9, SIR W (1) >9,..., S\R n (L) > 9) 
P(SIR > 9) 



1 



ES(-l) £ (T)exp 4 ( 




l-Pb 



where the expression in the numerator follows from [18 1. Moreover, for the special case of 
L = 1, P(SIRi > 6>|SIR >9) = {1- P b ) exp (2X9 2 / a d 2 {l - P b ) W^csc (f )) EH. ■ 

Hence using P out — P h + (1 - P,)Pf a ii|no backoff, we get the transmission capacity C — A(l - 
P ou t)R for CSMA MAP by combining Theorem |5] and [6] Finding the closed form expression 
for Pfaii| no backoff derived in Theorem |6] is quite challenging. An upper bound on the Pf ai i| no backoff, 
however, can be found using the FKG inequality EUll as follows. 

Definition 3: Let J 7 , V) be the probability space. Let A 6 J 7 , and 1^ be the indicator 
function of A. Event AG Jis called increasing if 1a (w) < 1,4(0/), whenever w < u/ for some 
partial ordering on to. The event A is called decreasing if its complement A c is increasing. 

Lemma 4: (FKG Inequality ) If both A, B E J 7 axe increasing or decreasing events then 
P{AB) > P(A)P(B) COD . 

Lemma 5: The outage probability of CSMA MAP P out < 1 - (1 - P b ) L+1 . 
Proof: Since SI R„(t) is decreasing function of the number of interferers, by considering cu = 
(ai, d2, . . • , ) where for m E N, a m = 1 if transmitter T m is active, and otherwise, it follows 
that the success event {SIR n (t) > 9} is a decreasing event. Hence, from the FKG inequality, 
P(SIR n (0) > 9, SIR n (l) > 9, . . . , SIR n (L) > 9) > P(SIR > 9) L+1 , since SIR n (t) is identically 
distributed for any t. Hence, P m \ no backof f < 1 - (1 - P b ) L , and P out < 1 - (1 - P b ) L+1 . ■ 

Discussion: In this section, we considered the CSMA MAP for an ad hoc network with energy 
harvesting nodes. We derived expressions for back-off and outage probability for the CSMA 
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MAP, thereby characterizing the transmission capacity. We showed that depending on the rate of 
energy arrival p, back-off probability can be written in closed form or can be expressed in terms 
of Lambert's function. We also derived an exact expression (and a simplified lower bound) for 
the outage probability, to characterize the transmission capacity with CSMA MAP. 

V. Simulations 

In this section, we present some numerical examples to illustrate our theoretical results. In all 
simulations we use energy arrival rate p = .5, a = 3, 9 = 2, and d = 2, such that A max = .023, 
except for Fig. [5] where a = 3, 9 = 1, and d = 1, and \ max = .2632 is used. In Figs. [4] and [51 we 
plot the transmission capacity for ALOHA MAP for B = oo with > p and <p = 0.5, 
respectively, from which we can see that for < p, the optimal transmission probability is 
q* = ^aai, while in the other case p < q* < 1, as derived in Theorem [T] In Figs. [6] and [7J 
we plot the transmission capacity for ALOHA MAP with B = 1, and B = 5, respectively. 
From Figs. SI |6] and [71 we can see that as B increases q* goes from 0.3 for B = 1 to 0.23 
for B = oo for fixed set of parameters. Finally, in Fig. [8l we plot the back-off probability with 
CSMA MAP as a function of A. Following Theorem [51 we see that for A = .01 and .035 for 
which ~ ]n ( p ^ Xmax > p 5 back off probability is equal to 1 — exp~( x ^), while for A = .05 and 

-If U / r(l-P b )A -\ 

.1, where — n w rnax < ^ satisfies the equation P b = 1 — exp v Xma * > . 

VI. Conclusions 

In this paper we considered ALOHA and CSMA MAP for an ad hoc network, and derived op- 
timal transmission probability for ALOHA MAP, and back-off and outage probability expressions 
for CSMA MAP, when each node in the network harvest energy from nature. We characterized 
the dependence of system throughput on the energy arrival rate, and derived system parameters 
for optimal performance. In this work, we assumed that each transmitter attempts to transmit with 
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same probability irrespective of the current energy state. With finite battery capacity, it makes 
more sense to transmit aggressively in high energy states and vice-versa. Analyzing energy aware 
transmission strategies remains an important problem to solve in future with energy harvesting 
nodes. 

Appendix A 

Proof of Lemma [T] We first consider the case of B = 1. For B — 1, fi(q) = — 22 — . Hence 

2 

/iw) = (p+q-pq) 2 an( * hence f[(q) > for q £ [0, 1]. For B > 1, for p ^ q, we next show that 
f' B (q) > for q £ [0,p) U (p, 1]. Let c := • F° r B > 1, g ^ p we compute the first 

derivative of f B (q) as f B (q) = > for q £ [0,p) U (p, 1]. 

In the interest of space we do not provide intermediate steps and only write the final answer. 
Also note that for q = p with B > 1, f B (q) = rq = B +f_ ■ It can be checked that f B (p — 5) < 
}b{jp) < /b(p + 8), for small 5 > 0. Hence is an increasing function of q G [0, 1]. 
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s = l-p(l-q)-q(l-p) 

Fig. 1. Transition state probabilities of infinite state birth-death Markov process. 
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s = l-p(l-q)-q(l-p),t=l-q(l-p) 

Fig. 2. Transition state probabilities of finite state birth-death Markov process. 



Distribution of packet arrivals to spatial nodes 
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Fig. 3. Packet arrival model for CSMA MAP. 
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Fig. 6. Plot of transmission capacity for B — 1 with p = 0.5. 
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Transmission Capacity for B=5 vs q for^.=0.1 , X = .02327 and p=.5 
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Fig. 7. Plot of transmission capacity for B — 5 with p = 0.5. 




Fig. 8. Backoff probability as a function of A. 



